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The  purpose  of  this  paper  is  to  develop  control  theoretic  concepts  for  distributed  control 
of  gas  turbine  engines,  and  develop  a  dynamic  engine  model  incorporating  distributed 
components  in  compressor  dynamics,  engine  cycles,  and  engine  control.  The  latest  results  in 
distributed  control  combined  with  adaptive  control  theory  are  extended  for  turbofan  engine 
distributed  control.  Concepts  and  architectures  for  distributed  control  are  developed  that 
create  tangible  benefits  from  the  distribution  of  closed-loop  feedback  around  the  engine. 


I.  Introduction 

Engineers  and  researchers1  agree  that  the  next-generation  engine  controls  need  to  move  to  a  distributed 
architecture  in  order  to  increase  flexibility  through  common  standards,  improve  redundancy  properties  by 
improving  the  overall  system  topology,  and  enable  component  self-diagnosing  and  other  benefits  of  ’smart’ 
actuators  and  sensors,  such  as  reduced  harness  weight.  In  this  vision,  distributed  computing  in  the  smart 
components  allows  for  localization  of  A/D  conversion  and  signal  processing,  supports  open  standards  and 
modularity,  and  provides  an  opportunity  for  self-diagnosis.  This  however  is  still  a  somewhat  limited  use  of 
distributed  computing  capabilities:  The  smart  components  neither  close  any  local  control  loops,  nor  perform 
any  functions  effecting  the  stability  or  performance  of  the  engine.  Thus  the  Full  Authority  Digital  Engine 
Controller  (FADEC)  still  remains  the  central  arbiter  of  the  engine’s  dynamic  behavior,  performance,  and 
reliability. 

The  problem  we  will  address  in  this  effort  is  the  need  to  go  beyond  the  limited  application  of  distributed 
computing.  This  need  arises  because  the  full  potential  of  a  distributed  architecture  cannot  be  attained 
unless  the  control  algorithms  themselves  are  distributed.  For  instance,  if  the  control  laws  are  not  distributed 
the  dependence  on  the  FADEC  remains  high,  and  system  reliability  can  only  be  insured  through  many 
redundant  components  and  interconnections.  Information- flow  and  redundancy  requirements  are  still  based 
on  a  centralized  controller,  and  the  potential  for  complex  and  hard-to-modify  centralized  code  remains. 
Furthermore,  the  benefits  of  adaptation,  robustness,  and  self-repair  at  the  component  level  are  not  attained. 

*PhD  Candidate,  School  of  Aerospace  Engineering,  Mehrdad.Pakmehr@gatech.edu,  AIAA  Student  Member. 

^MS  Student,  School  of  Industrial  and  Systems  Engineering,  marion.mounier@gatech.edu 

-I- Propulsion  Engineer,  Aurora  Flight  Sciences  Corporation,  nfitzgerald@aurora.aero. 

§  Propulsion  Engineer,  Aurora  Flight  Sciences  Corporation,  gkiwada@aurora.aero. 

^Principal  Investigator,  Aurora  Flight  Sciences  Corporation,  jpaduano@aurora.aero,  AIAA  Member. 

II  Professor,  School  of  Aerospace  Engineering,  feron@gatech.edu,  Associate  Fellow  AIAA. 

** Senior  Aerospace  Engineer,  Air  Force  Research  Laboratory,  al.behbahani@wpafb.af.mil,  Associate  Fellow  AIAA 

Copyright  ©  2009  by  the  American  Institute  of  Aeronautics  and  Astronautics,  Inc.  The  U.S.  Government  has  a  royalty-free 
license  to  exercise  all  rights  under  the  copyright  claimed  herein  for  Governmental  purposes.  All  other  rights  are  reserved  by  the 
copyright  owner. 


1 


American  Institute  of  Aeronautics  and  Astronautics  Paper  AIAA-2009-5532 


Such  properties  are  best  achieved  through  feedback  control  at  the  component  level  and  have  the  potential 
to  improve  overall  system  reliability. 

Previous  work  on  distributing  the  control  algorithms  in  gas  turbine  engines  has  not  taken  advantage  of 
recent  progress  in  adaptive  control  algorithms.  An  adaptive  controller  requires  little  or  no  a  priori  information 
about  the  unknown  parameters,  improves  its  performance  as  it  adapts,  and  is  likely  to  use  lower  gains  since 
adaptation  extends  on  to  the  level  that  is  necessary  instead  of  to  a  conservative  pre-determined  level.  The 
probable  advantages  of  adaptive  control  are  particularly  applicable  to  distributed  control  in  gas  turbine 
engines  as  they  can  probably  help  overcome  communication  imperfections  such  as  lag  or  packet  drop  inherent 
in  a  networked  control  architecture. 

To  determine  the  feasibility  of  a  distributed  adaptive  engine  control  approach,  a  systematic  evaluation  is 
needed  of  the  stability  and  performance  characteristics  that  would  result  if  not  only  sensors  and  actuators 
but  entire  subsystems  of  the  engine  become  ’smart’;  that  is,  if  distributed  computing  is  used  at  the  local 
level  and  only  coordinated  by  the  FADEC.  Such  an  architecture  must  be  studied  in  the  context  of  noisy, 
band-limited  and  delayed  communications  between  the  subsystem  controllers.  Furthermore,  the  impact  of 
varying  operating  conditions  on  the  performance  at  the  global  level  of  the  hierarchical  control  system  must 
be  addressed.  Finally,  the  potential  for  instabilities  due  to  the  interaction  of  separate  controllers  when  delays 
and  synchronization  issues  arise  must  also  be  addressed. 

Distributed  control  is  a  broad  term  that  encompasses  many  different  levels  of  control  delegation  from  the 
traditional  FADEC.  As  the  signal  processing  and  control  laws  are  distributed  further  there  are  also  different 
configurations  in  which  engine  components  can  be  grouped  for  control.  Some  sample  distributed  control 
architectures  are  compared  with  a  FADEC  in  figure  1.  Benefits  from  distributed  control,  such  as  sensor 
modularity,  weight  reduction,  and  life-cycle  cost  reduction  have  been  discussed  in  previous  works.1  3  An 
appealing  configuration  is  the  Partially  Distributed  control  architecture,  where  there  are  local  controllers  with 
some  authority  but  the  entire  engine  system  is  still  governed  by  a  central  supervisory  controller.  Partially 
distributed  control  has  many  of  the  benefits  of  fully  distributed  control  but  retains  a  central  supervisor 
to  communicate  with  the  operator  and  handle  some  system  level  tasks  such  as  engine  start  up.  Partially 
distributed  control  architecture  enables  a  new  engine  development  paradigm  infeasible  with  less  distributed 
control  schemes. 


Figure  1.  Control  Architecture  Diagrams  for  a  Gas  Turbine  Engine 
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The  new  engine  development  scenario  envisioned  here  has  gas  turbine  engine  cores  being  utilized  in  a 
more  modular  manner  similar  to  the  current  use  of  internal  combustion  (IC)  engines.  As  an  illustrative 
example  consider  a  simple  IC  engine  Helicopter.  The  IC  engine  is  purchased  with  its  own  governor  or  Engine 
Control  Unit  (ECU)  and  connected  via  a  transmission  to  the  rotor,  a  variable  load.  The  operator  has  control 
over  the  engine  RPM  as  well  as  the  engine  load  (via  the  rotor).  Applying  this  analogy  to  a  gas  turbine 
instead  of  a  helicopter  the  IC  engine  represents  a  gas  generator  core  (High  Pressure  Compressor  (HPC), 
combustor,  High  Pressure  Turbine  (HPT)),  the  rotor  load  represents  either  a  shaft  driven  device  or  even 
a  new  spool  (Fan,  Low  Pressure  Turbine  (LPT)),  and  the  operator  is  a  supervisory  controller.  Creating 
a  separate  engine  controller  for  the  gas  generator  has  strong  commercial  applications,  both  in  large  scale 
commercial  gas  turbine  design  and  small  scale  UAV  development. 

Commercial  manufacturers  of  Gas  Turbine  engines  rarely  design  all  new  engine  centerlines,4  the  lifespan 
of  successful  engine  families  lasts  decades.  Many  of  the  new  engines  designed  in  a  family  are  based  on  an 
existing  engine  core,  primarily  due  to  cost  and  reliability  concerns.  The  high  pressure  compressor  and  turbine 
contain  the  highest  performance,  and  therefore  most  expensive,  components.  Engine  core  designs  may  move 
from  military  turbojets  into  commercial  turbo- fans  and  turbo-props.5  In  the  case  of  the  extremely  popular 
CFM56  (figure  2)  a  GE  F101  engine  was  used  as  the  core  for  a  different  company’s  engine. 


i - VS  Sneema 

tan  '^55-  *  SAFftAN  Group 


Figure  2.  CFM56  engine  developers  by  component6 

Distributed  adaptive  engine  control  could  enable  ’’plug  and  play”  development  of  entire  families  of  engines. 
In  the  distributed  engine  control  vision  engine  cores  could  be  purchased  with  onboard  subordinate  controllers 
ready  for  integration  into  a  larger  engine,  whereas  the  CFM56  FADEC  was  developed  independently  for  the 
integrated  engine.  Structuring  engine  control  in  such  a  distributed  fashion  would  increase  compatibility 
between  different  engine  manufacturers  and  reduce  development  time  and  cost  for  new  engines. 

A  similar  niche  is  occurring  in  UAV  development,  where  small  gas  turbines  are  being  used  to  power 
a  variety  of  different  lift/thrust  devices.  UAV  development  programs  rarely  have  the  resources  for  serious 
engine  redevelopment  and  therefore  must  select  from  a  limited  number  of  commercial  off  the  shelf  (COTS) 
engines.  In  the  case  of  small  gas  turbines  these  COTS  engines  are  generally  designed  for  missile-turbojet  or 
power  generator  applications,  while  the  UAV  designer  may  want  to  use  the  engine  core  in  a  turbo-prop  or 
turbofan  application.  Successful  development  of  adaptive  distributed  control  for  this  class  of  engines  would 
allow  UAV  designers  to  purchase  engines  with  onboard  controllers  and  mate  them  with  their  own  proprietary 
fan/prop  sections  without  having  to  design  a  new  control  system  from  scratch. 

In  this  paper  we  demonstrate  the  feasibility  of  a  partially  distributed  control  scheme  with  separate 
controllers  on  the  engine  core  and  fan,  where  the  controllers  are  linked  by  a  supervisory  controller.  This 
scheme  is  representative  of  the  situation  encountered  in  VTOL  UAV  design  and  the  design  of  new  turbo¬ 
props  and  variable  pitch  turbofans  by  the  large  commercial  gas  turbine  manufacturers.  For  future  we  will 
develop  the  partially  distributed  controller  further  to  cover  safe  performance  during  non  standard  operations 
(including  sensor  failure  etc.),  culminating  in  a  static  engine  test  of  a  small  turbo-prop  engine  running  the 
developed  distributed  adaptive  controller. 
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II.  Engine  Model 


A.  Overview 

For  the  purpose  of  the  distributed  control  analysis  we  use  a  simplified  dynamic  model  of  a  generic  turboshaft 
engine  driving  a  lift  fan,  depicted  in  Figure  3.  Fuel  is  provided  to  the  core  of  the  engine,  which  is  comprised 
of  a  compressor,  combustor,  and  turbine.  The  turbomachinary  components  are  connected  by  a  shaft  that 
also  provides  torque  to  the  fan  by  way  of  a  reduction  gearbox.  The  thrust  from  the  fan  can  be  controlled  by 
varying  the  pitch  angle  of  the  fan  blades,  or  by  varying  the  exit  area  of  fan  duct. 

i 


Figure  3.  Turbofan  engine  diagram 

The  dynamics  of  thrust  and  efficiency  are  modeled  using  the  lumped-parameter,  first-principles  approach 
common  to  standard  gas  turbine  textbooks.7,8  The  performance  of  the  core  will  be  represented  by  the 
power  of  the  turbomachinary  components,  which  can  be  described  as  functions  of  the  inlet  temperatures 
and  pressures  through  the  engine,  the  air  flow  through  the  core,  pressure  ratios  across  the  components, 
component  efficiencies,  and  temperature  in  the  combustor.  The  fan  performance  can  be  determined  from 
the  airflow  through  the  fan,  the  fan  pressure  ratio,  and  inlet  temperature  and  pressure. 

The  model  of  the  engine  system  for  this  phase  uses  the  following  assumptions  and  simplifications: 

•  Only  the  dynamics  of  the  spool  will  be  considered.  Other  low  speed  dynamics,  like  heat  transfer  from 
the  gas  path  to  metal,  or  higher  speed  dynamics,  such  as  acoustics,  volume  dynamics,  and  combustor 
heat  release  dynamics,  are  ignored  in  the  current  analysis.  The  essential  role  of  the  engine  controller 
can  be  assessed  from  this  perspective,  and  in  future  phases  of  the  project  more  detailed  dynamics  can 
be  modeled  if  necessary. 

•  The  analysis  will  assume  sea  level  static  operation  of  the  engine  system,  typical  of  a  hovering  condition 
for  a  VTOL  UAV. 

•  The  turbine  expands  the  exhaust  gases  of  the  core  perfectly  to  ambient  pressure,  such  that  all  possible 
work  is  extracted  from  the  flow.  The  exhaust  of  the  core  is  assumed  to  have  no  impact  on  the  thrust 
produced  by  the  fan. 

•  Non-ideal  efficiencies  are  assumed  for  the  fan,  compressor,  combustor,  and  turbine,  but  all  other 
components  are  assumed  to  operate  ideally.  In  particular,  there  is  no  pressure  loss  in  the  inlets  the 
core  and  fan,  and  the  fan  nozzle  is  assumed  to  allow  perfect  expansion  of  the  fan  airflow  back  to 
ambient  pressure. 

•  The  core  is  modeled  as  a  simple  gas  generator  with  no  secondary  or  bleed  flows  for  turbine  cooling  or 
other  uses. 

•  The  mass  and  temperature  of  the  fuel  is  ignored  in  the  calculations.  The  fuel  is  simply  modeled  as 
adding  heat  to  the  thermodynamic  cycle  in  accordance  with  its  mass  flow  rate  and  heat  of  combustion. 

•  The  calculations  assume  that  the  specific  heat  of  the  air  going  through  both  the  core  and  the  fan  is 
independent  of  temperature. 
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B.  Detailed  Description 

In  the  model  presented  here,  the  performance  of  the  engine  system  is  a  function  of  several  lumped-parameter 
thermodynamic  variables.  For  the  core,  these  are  the  corrected  mass  flow  through  the  compressor  (Wc), 
the  compressor  pressure  ratio  (7rc),  and  the  ratio  of  turbine  inlet  stagnation  temperature  to  the  compressor 
inlet  stagnation  temperature  (T4/T2).  For  the  fan,  the  key  thermodynamic  parameters  are  the  corrected 
mass  flow  (Wf)  and  the  pressure  ratio  (7 if).  For  a  given  inlet  stagnation  temperature  and  pressure  (T2  and 
P2),  the  thrust  and  the  fuel  consumption  of  the  engine  can  be  set.  The  spool  speed,  TV,  is  the  mechanical 
constraint  that  links  the  thermodynamic  condition  of  the  individual  components  to  each  other. 

The  dynamics  of  the  spool  can  be  described  from  Newton’s  second  law  by  the  sum  of  the  torques  produced 
by  the  turbomachinary  components  (T),  as  well  as  the  overall  rotational  inertia  of  the  system  (/). 

Iu  =  Tt-Tc-Tf .  (1) 

Here,  u  is  expressed  in  rad/s  and  the  subscripts  t,  c,  and  /  represent  the  turbine,  compressor,  and  fan 
respectively.  Noting  that  the  power  output  of  each  of  those  components  is  defined  by  V  =  Tuj  and  expressing 
the  key  parameters  as  a  fraction  of  the  design  value,  Equation  1  can  be  rewritten  as 

_  f  60  Vdes  \  /  Vt  Vt,des  Vc  Vc,des  Vf  \  ,  (  N  \  ,^\ 

Ndes  ~  WlNlJ  \  V t,des  V des  V c,des  V 'des  Vf^des)  \  ^des  ) 

where  the  rotational  speed  (N)  has  been  expressed  in  rpm  instead  of  rad/s.  The  expression  of  the  performance 
parameters  as  a  ratio  of  their  design  values  allows  the  analysis  to  circumvent  the  knowledge  of  several  detailed 
parameters  in  the  cycle.  For  instance,  the  gear  ratio  of  the  gearbox  is  not  required  in  this  representation. 
The  power  output  of  each  of  those  components  is  determined  by  the  thermodynamic  matching  of  the  engine. 

The  power  requirement  of  the  compressor  is  equal  to  the  enthalpy  increase  from  station  2  to  station  3  in 
Figure  3,  which  can  be  expressed  as  a  function  of  mass  flow,  pressure  ratio,  adiabatic  efficiency,  and  inlet 
conditions.  Written  as  a  fraction  of  the  design  power, 


Vc  =  W2  P2  I  T2  J)cAes  7TC  7  -  1 

Vc^des  W2,des  P2,std\  ^2,std  Vc  J^T  _i’ 
V  nc,des  1 


(3) 


assuming  constant  ratio  of  specific  heats,  7  =  cp/cv.  The  mass  flow  is  expressed  as  corrected  flow,  or 
equivalent  flow  at  standard  conditions  W  =  rriyjT /Tstd/ (P / Pstd)-  P  is  stagnation  pressure  and  77  is  adiabatic 
efficiency. 

The  fan  power  consumption  is  exactly  analogous  to  that  of  the  compressor, 


Vf 

wf 

P2 

/  T2 

Vf,des 

7-l 

v  -1 

Vf^des 

^Vf,des 

P2,std  ^ 

1  T2,std 

Vf 

7~1 

nLdes  ~  1 

The  power  provided  by  the  turbine  is  equal  to  the  energy  extracted  as  the  gas  is  expanded  from  the 
combustor  exit  temperature  and  pressure  back  to  ambient  conditions.  Assuming  static  conditions,  the 
pressure  ratio  across  the  turbine  is  therefore  set  by  the  compressor  pressure  ratio  and  the  burner  pressure 
loss  (71771777  =  1).  Assuming  that  the  mass  flow  through  the  turbine  is  equal  to  that  of  the  compressor,  the 
power  output  can  be  calculated  as  a  function  of  the  turbine  inlet  temperature. 


Vt 

V t,des 

The  turbine  inlet  temperature  is  determined  by  the  amount  of  fuel  entering  the  combustor.  Ignoring  the 
mass  of  the  fuel  and  the  fuel  temperature,  the  inlet  stagnation  temperature  to  the  turbine  can  be  calculated 
from  a  simple  energy  balance. 

rhfhf  =  rhccp  (T4  -  T3) .  (6) 

As  written,  this  equation  assumes  that  all  of  the  energy  of  the  fuel  is  translated  to  heat  in  the  combustor 
regardless  of  equivalence  ratio  or  adiabatic  flame  temperature  of  the  fuel  air  reaction.  In  the  actual  model 


Wc  P2  /  T2  rjt  1  (7777) 


-(7-1) 


w, 


c,des  2 ,std  y  -^2 ,std  Vt,des  ^  (ltc  des^b  des) 


-(7-1)  ' 


(5) 
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code,  an  upper  limit  is  placed  on  T4  near  the  max  adiabatic  flame  temperature  to  prevent  the  simulation 
from  adding  a  limitless  amount  of  energy  into  the  system.  Expressed  as  a  ratio  of  design  values,  the  effect 
of  fuel  flow  on  turbine  inlet  temperature  can  be  expressed  as 

F1 fuel  _  Wc _ t2  Tc  ^ 

FlfUel,des  Wc,des  3,des  —  Tc  des 

-L2,des  ’ 

Note  that  the  temperature  ratio  across  the  compressor  (rc)  is  known  from  the  compressor  pressure  ratio  (7 rc) 
and  adiabatic  efficiency  (rjc). 

The  power  expressions  from  Equations  3,  4,  and  5  are  linked  together  and  to  the  spool  speed  through 
mass  flow  conservation,  as  well  as  by  the  individual  component  performance  maps.  In  the  core,  the  mass 
flow  is  set  by  the  throttling  effect  of  the  turbine,  where  the  nondimensional  flow  through  the  component  is 
set  by  the  pressure  ratio  across  it.  The  inlet  conditions  of  the  turbine  can  be  expressed  as  functions  of  the 
compressor  and  combustor  performance  as  follows: 


rhc^/RT^ 

P4A4 


=  FP4^t)  = 


m, 


VBT2 


P0A0  Fs  E± 

J  22  p2  p3 


(8) 


For  a  fixed  design,  compressor  and  turbine  inlet  areas  (A2  and  A4)  are  considered  constant,  so  they  fall 
out  when  the  equation  is  expressed  as  a  ratio  of  design  values. 


7TC  7 Tb  FP4  /  T^dea  =  Wc 

^ c,des  'K b,des  FP4 ,des  y  ^2 ,des  ^Fc,des 

To  simplify  the  calculation,  the  model  assumes  that  the  turbine  is  thermodynamically  choked,  so  FP4  is 
constant  and  equal  to  the  design  value.  This  is  a  typical  assumption  made  in  turbine  engine  models,  and  is 
usually  invalid  only  at  engine  settings  near  idle. 

For  the  fan,  the  feature  that  controls  the  mass  flow  is  the  nozzle  exit  area,  Ae.  Since  the  nozzle  is  not 
choked,  the  expression  for  the  conservation  of  mass  flow  must  take  the  flow  parameter  variation  into  account. 
Writing  fan  nozzle  mass  flow  as  a  function  of  fan  inlet  and  operating  conditions,  we  have 


rhf  \JRTe  _  /  \  _  Fif  VRT2  yfc  1 

p  A  ~  r  ^e\Pf)  —  p  A  pp  a  • 

FeAe  P2A2  ^  ^ 

Dividing  the  parameters  in  this  equation  by  their  design  values,  we  have 

7 r/  Ae  FPe{nf )  =  Wf  I  rf 

'Kfides  Aejdes  FPe(jT f,des)  VVf^des  \  7 ~f,des 


(10) 


(11) 


The  nozzle  flow  parameter  characteristic,  FPe(7Tf),  can  be  calculated  for  static  conditions  assuming  ideal 
expansion  to  ambient  pressure  as 


FP  = 


(12) 


Component  maps  provide  the  final  link  from  thermodynamic  performance  to  mechanical  state.  For  the 
turbine,  we  assume  that  the  adiabatic  efficiency  is  constant  and  equal  the  the  design  value  for  the  range  of 
conditions  to  be  simulated.  For  the  fan  and  compressor,  the  pressure  ratio  and  efficiency  of  each  component 
is  expressed  as  a  function  of  the  spool  speed  and  airflow,  both  corrected  to  inlet  conditions.  The  change  in 
incidence  angle  of  the  fan  blades  from  the  design  condition (Ai)  is  also  input  to  the  fan  map. 


TTc  Vc 
c,des  Vc,des 


Compressor  Map 


(13) 


*7  Vf 

_^f,des  Vf,des_ 


FanMap 


wf  n  rw  A.\ 

Wftdes' NdJ\l  Tstd  / 


(14) 
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For  the  purposes  of  this  Phase  of  investigation,  generic  maps  have  been  generated  for  the  fan  and  compressor 
based  on  assumed  0/t/?  characteristics  and  other  aspects  of  compressor  theory.  As  the  model  is  improved, 
these  generic  maps  can  be  replaced  with  ones  representing  the  performance  of  specific  compressors  and  fans. 

For  a  given  spool  speed,  the  seven  unknown  parameters  of  the  system  (Wc,  Wf,7rc,7Tf,r]c,rif,  and  T4/T2) 
can  be  calculated  from  a  system  of  seven  equations,  Equations  7,  9,  11,  13  and  14.  Note  that  13  and  14 
represent  two  equations  each.  The  solution  of  that  system  of  equations  provides  all  the  necessary  information 
to  determine  the  dynamics  of  the  system  as  expressed  in  Equation  2. 

Lastly,  the  performance  output  of  the  engine  system,  namely  thrust  and  fuel  consumption,  can  be  calcu¬ 
lated  from  the  solved  system  parameters.  At  static  conditions,  the  thrust  from  the  fan  can  be  approximated 
from  first  principles  as  F  =  rrifUe ,  where  ue  is  the  velocity  of  the  air  exiting  the  fan, 


F 

AePT2 


=  7  M2e 


As  a  function  of  design  parameters,  this  becomes 


27 

7-1 


7—1 


i). 


F  _  Ae  P2  V  -1 

zp  A  p  7-1 

r  des  ^-e^des  Jr2,des  7  1 

7TJ‘,des  1 

Specific  fuel  consumption  is  follow  as  the  ratio  of  fuel  flow  to  thrust 


SFC  m  t  mfuel  \  (  F  \ 

P  FF 'des  \Alfuel  ,des  /  \Pdes  J 


(15) 


(16) 


(17) 


C.  Linearized  Engine  Model 

Nonlinear  engine  model  can  be  represented  as 

x(t)  =  f(x(i),u1(t),u2(i),u3(t)),  /lgx 

y(t)  =  g(x(t),u2(t),u3(t)), 

where  x(t)  is  the  spool  speed  and  the  only  state  of  the  system,  u\(t)  is  the  fuel  flow  control  input,  u2(t) 
is  the  fan  exit  area  control  input  and  u3(t)  is  the  fan  vane  angle  control  input.  Control  inputs  ranges  are 
ui(t)  G  [0.3,  1.2],  u2(t)  G  [0.8,  1.2]  and  u3(t)  G  [—10,  10] (deg)  which  are  imposed  by  physical  limitations. 
Linearizing  this  model  around  design  point  ( Xdes  =  1  ,Udes  =  [1,  1,  0]T)  we  have 

5x(t)  =  a.Sx(t )  +  bi.Sui(t)  +  62A  ^2  (t)  +  bs.Su^t),  (19) 

which  also  can  be  written  as 

Sx(t)  =  a.Sx(t )  +  b.Su(t ),  (20) 

where  Sx  =  x  —  5u  =  u  —  Udes-  For  the  purpose  of  closing  the  loop  in  the  fan  subsystem,  we  defined 

the  thrust  as  the  output  and  assumed  we  can  estimate  it  as  a  linear  function  of  the  state,  fan  vane  angle, 
and  fan  exit  area  control  inputs. 

5y(t)  =  c.5x(t )  +  d\.5u\(t)  +  d2.5u2(t)  +  ds.Su3{t ),  (21) 

where  Sy  =  y  —  ydes ?  Su  =  u  —  Udes ?  where  the  design  value  for  output  is  ydes  =  1-  Constants  are  a  = 
—4.63,  bi  =  1.78,  b2  =  0.31,  b3  =  2.96,  c  =  1.921,  d\  =0,  d2  =  0.215,  d3  =  —1.44.  These  values  has  been 
found  by  numerical  linearization  of  the  engine  dynamics. 


III.  Distributed  Linear  Control 

For  distributed  control  purpose,  we  divided  the  model  into  two  subsystems  which  are  the  core  engine  and 
the  fan.  Distributed  turbofan  engine  model  is  represented  in  figure  4. 2 

The  highlighted  subsystems  are  candidates  for  simplified,  distributed  control  structure.  In  our  engine 
model,  system  I  is  analogous  to  the  ’’Inlet  Fan”,  and  system  II  is  analogous  to  the  core  engine  (i.e.  Com¬ 
pressor,  Combustor,  Turbine  Nozzle). 
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Figure  4.  Schematic  of  turbofan  engine  subsystems2 


In  this  section  we  design  linear  controllers  for  two  subsystems  based  on  classical  control.  The  controllers 
will  be  I  and  PI  controllers.  The  figures  of  merit  which  we  use  to  quantify  the  performance  are  percent 
overshoot  (PO)  and  settling  time  (ts)  for  the  transient  response  and  steady  state  error  (ess)  for  the  steady 
state  response  of  the  system. 

A.  Core  Controller  Design  Strategy 

Here,  we  need  x(t)  track  a  given  smooth  trajectory  rc(£),  while  all  other  signals  remain  bounded. 

introducing  the  integral  error  of  spool  speed  as  an  additional  state  into  the  system  dynamics,  we  have 

x(t)  =  5xi(t)  —  rc(t),  x(0)  =  0.  (22) 

The  new  state  will  be  ^ 

5x2 (t)  =  xj(t)  =  /  x(r)dr.  (23) 

Jo 

The  model  which  we  use  to  control  the  engine  spool  speed  is  developed  assuming  5u2(t)  =  5us(t)  =  0 

5±i(t)  =  a.5xi(t)  +  bi.5ui(t),  &tl(0)  =  Sx2(  0)  =  0,  ^4) 

Sx2(t)  =  Sx!(t)  -  rc, 

where  5x\{t)  is  the  spool  speed  and  Sx2(t)  is  the  integration  of  spool  speed  tracking  error. 

The  linear  control  logic  for  fuel  flow  is 

5ui{t)  =  kpf.6xi(t)  +  kif.5x2{t).  (25) 

B.  Decentralized  Linear  Controller  Design  Strategy 

Decentralized  control  structure  we  present  here  has  two  subsystems  (i.e.  core  engine  and  the  fan)  and  we 
design  two  linear  controllers  for  each  subsystem.  Fuel  flow  is  the  control  input  for  the  core  subsystem  and 
vane  angle  is  the  control  input  for  the  fan  subsystem. 
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1.  Subsystem  I:  Fan  Control  Design 

Here,  we  need  y(t)  to  track  a  given  smooth  trajectory  yc(t),  while  all  other  signals  remain  bounded. 

Introducing  the  integral  error  of  the  output  as  an  additional  state  into  the  system  dynamics,  we  have 


where 


y(t)  =  Sy(t)  -yc(t), 


Sy(t)  =  c.5xi(t)  +  d3.Su3(t), 


is  the  linear  estimation  of  the  turbofan  thrust. 

The  new  state  is 

Sx3(t)  =  yi{t )  =  f  y(r)di 
Jo 

Hence  the  control  logic  for  fan  vane  angle  is 

^3  (t)  =  kiaSx3(t ), 


(26) 

(27) 


(28) 


(29) 


where  Sx3(t)  is  the  integration  of  output  tracking  error.  Note  that  d3  <  0,  hence  we  need  a  negative  gain. 


2.  Subsystem  II:  Core  Engine  Control  Design 

The  control  for  this  subsystem  is  similar  to  what  we  explained  in  the  ’core  control  design  strategy’  section. 
We  construct  the  rc(t)  signal  from  yc(t )  and  Su3(t)  using  the  thrust  estimation  equation.  The  following 
relation  is  used  to  reconstruct  rc(t)  from  yc(t ) 

rc(t)  =  [: yc{t )  -  ydes  -  d3(u3  -  u3des) (tt / 180)\ / c  +  xdes.  (30) 

The  decentralized  control  structure  in  which  we  construct  the  rc(t)  signal  from  yc(t)  and  Su3(t)  using  the 
thrust  estimation  equation,  is  shown  in  figure  5.  The  signals  are  shown  with  solid  lines  and  the  mechanical 
interconnection  is  shown  with  a  dashed  line. 


Figure  5.  Decentralized  linear  control  structure,  rc{t)  signal  depends  on  yc(t) 


C.  Core  Engine  Linear  Control  Results 

In  this  simulation  we  assumed  that  torque  is  a  function  of  time  and  defined  it  as  T(t)  =  1  +  0.1  sin (1.3(1  + 
sin(4t))t).  The  linearized  nominal  system  for  this  simulation  is  x(t)  =  a.Sx(t)  +  b.5ui(t)  =  —  1 .8Sx(t)  + 
1.7460^1  (£). The  simulation  results  for  linear  control  of  core  engine  are  presented  here. 
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The  control  parameters  and  initial  conditions  for  all  of  these  cases  are 


kpf  9 1  kif  75, 
x(0)  =  1.0,  T(0)  =  1.0. 


In  this  simulation  rc(t )  is  a  step  signal  which  defines  the  reference  spool  speed  to  change  from  1  to  0.95. 
Simulation  results  are  shown  in  figures  6  to  9.  Figure  6  shows  spool  speed  history  for  nominal  and  reference 
plant  and  also  state  tracking  error  history.  It  can  be  seen  that  state  follows  the  reference  trajectory  very 
closely.  Figure  7  shows  the  history  of  the  load  on  the  core  engine.  Figure  8  shows  fuel  flow  history  as  control 
input.  Figure  9  shows  excess  power,  core  power,  SFC,  turbine  temperature,  core  pressure  ratio,  and  core 
airflow  time  histories. 


spool  speed 


Load  history 


Figure  6.  Spool  speed  and  integration  of  state  tracking 
error  histories 


Figure  7.  Load  history 


fuel  control  input  history 


Core  engine  power 


Turbine  Temperature 


time(sec) 

Core  engine  pressure  ratio 


time(sec) 

Core  engine  air  flow 


Figure  8.  Fuel  flow  control  input  history 


Figure  9.  Turbine  temperature,  SFC,  pressure  ratio, 
power  and  airflow  histories 


D.  Turbofan  System  Decentralized  Linear  Control  Results 

Here  we  show  the  results  of  the  simulation  for  the  decentralized  case  with  two  different  linear  controllers  to 
control  two  different  subsystems  of  the  turbofan  engine  (i.e.  engine  core  and  fan).  Fuel  flow  and  fan  vane 
angle  are  the  control  inputs. We  use  step  reference  signal  ( yc(t ))  for  the  output  thrust  to  track.  rc(t)  will  be 
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constructed  using  yc(t )  and  Sus{t).  The  decentralized  control  structure  that  we  use  in  this  section  is  shown 
in  figure  5. 

The  control  parameters  and  initial  conditions  for  all  of  these  cases  are 

kpf  —  9,  kif  =  75,  kia  =  10, 

x(0)  =  1.0.  (  ' 

In  this  simulation  yc(t )  is  a  step  signal  which  defines  the  reference  spool  speed  to  change  from  1  to  0.95. 
Simulation  results  are  shown  in  figures  10  to  13.  Figure  10  shows  the  history  of  the  load  on  core  engine,  and 
also  estimated,  actual  and  reference  thrust  time  histories.  Figure  11  shows  spool  speed  history  for  nominal 
and  reference  plant,  and  also  time  histories  of  integration  of  the  state  tracking  error  and  (i.e.  xi)  output 
tracking  error  (i.e.  yj).  The  state  follows  the  reference  trajectory  closely.  Figure  12  shows  fuel  flow,  vane 
angle  and  exit  area  histories  as  control  inputs.  Figure  13  shows  excess  power,  core  power,  SFC,  turbine 
temperature,  core  pressure  ratio,  fan  pressure  ratio,  core  airflow  and  fan  airflow  time  histories. 
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Figure  10.  Plant  actual  thrust,  estimated  thrust,  ref-  Figure  11.  Spool  speed,  and  integration  of  the  state 
erence  thrust,  and  load  histories  and  the  output  tracking  errors  histories 
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Figure  12.  Fuel  flow,  fan  exit  area  and  fan  vane  angle  Figure  13.  Power,  turbine  temperature,  SFC,  fan  and 
control  inputs  histories  core  pressure  ratio  and  airflow  histories 
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IV.  Distributed  Adaptive  Control 


Since  engine  dynamics  have  both  parametric  and  dynamic  uncertainties,  we  decided  to  use  adaptive 
control  technique  to  control  the  engine.  The  basic  idea  in  adaptive  control  is  to  estimate  the  unknown 
parameters  on-line  based  on  measured  system  signals,  and  use  the  estimated  parameters  in  the  control  input 
computations.  Since  adaptive  control  systems  are  inherently  nonlinear,  most  of  the  times,  their  design  and 
analysis  is  strongly  connected  with  Lyapunov  stability  theory. 

In  this  section  we  develop  a  distributed  (decentralized)  control  strategy  for  a  turbofan  engine  using  model 
reference  adaptive  control  (MRAC)  technique.  The  controller  design  methodology  is  mostly  based  on  the 
results  developed  in  references.9, 10 

A  model-reference  adaptive  control  (MRAC)  system  is  composed  of  four  parts:  a  system  containing 
unknown  parameters,  a  reference  model  for  specifying  the  desired  output  of  the  control  system,  a  feedback 
control  law  containing  adjustable  parameters,  and  an  adaptation  mechanism  for  updating  the  adjustable 
parameters. 

The  adaptation  mechanism  is  used  to  adjust  the  parameters  in  the  control  law.  The  objective  of  the 
adaptation  is  to  make  the  tracking  error  converge  to  zero.  The  difference  from  conventional  control  lies  in  the 
structure  of  this  mechanism.  The  main  issue  is  to  synthesize  an  adaptation  mechanism  which  will  guarantee 
that  the  control  system  remains  stable  and  the  tracking  error  converges  to  zero  as  the  parameters  are  varied. 
Adaptation  law  is  not  necessarily  uniquely  defined. 

Absence  of  stability  margin  metrics  and  lack  of  solid  theoretical  results  on  distributed  adaptive  systems 
with  communication  constraints  (i.e.  delay,  packet  drop,  etc)  are  the  key  areas  to  be  developed  in  distributed 
adaptive  control  systems. 

In  the  case  of  stability  margins  we  can  implement  the  methods  developed  by  C.  Cao  and  N.  Hov- 
akimyan11  14  which  present  adaptive  control  architecture  that  adapts  fast  and  ensures  uniformly  bounded 
transient  response  for  systems  both  signals,  input  and  output,  simultaneously.  This  architecture  has  a  low- 
pass  filter  in  the  feedback  loop  and  relies  on  the  small-gain  theorem  for  the  proof  of  asymptotic  stability.  The 
tools  from  these  papers  can  be  used  to  develop  a  theoretically  justified  verification  and  validation  framework 
for  adaptive  systems. 

In  case  of  the  distributed  adaptive  systems  with  delay  we  can  again  use  the  results  developed  by  C.  Cao 
and  N.  Hovakimyan15, 16  which  extend  the  results  from  their  previous  paper11  for  characterization  of  the 
time-delay  margin  of  closed  loop  systems. 

A.  Controller  Design  Strategy 

The  decentralized  adaptive  control  structure  for  the  turbofan  engine  model  is  shown  in  figure  14.  The  signals 
are  shown  with  solid  lines  and  the  mechanical  interconnection  is  shown  with  a  dashed  line. 


1.  Subsystem  II:  Core  Engine  Control  Design 


The  control  objective  for  subsystem  II  is  to  design  an  adaptive  controller  to  achieve  tracking  x(t)  — >  xm(t) 
as  t  — >  oo. 

Here  we  introduce  the  integral  error  in  spool  speed  as  an  additional  state  into  the  system  dynamics  and 
write 

x(t)  =  Sxi(t)  —  rc(t)f  x(0)  =  0.  (33) 


The  new  state  will  be 


6x2  (t)  =  xi(t)  =  /  x{r)dr. 

J  o 


(34) 


The  model  which  we  use  to  control  the  engine  spool  speed  is  developed  assuming  Su2{t)  =  Sus(t)  =  0 


6xi(t)  =  a.5xi(t)  +  bi.5ui(t),  fei(O)  =  fe2(0)  =  0, 

Sx2(t)  =  Sxi(t)  -  rc, 

where  6x\(t)  is  the  spool  speed  and  Sx2(t)  is  the  integration  of  spool  speed  tracking  error. 
Let  X(t)  =  [Sxi(t),  Sx2(t)]T .  The  system  dynamics  (51)  can  be  written  as 
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Figure  14. 


Decentralized  adaptive  control  structure  for  turbofan  engine 
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(36) 


Let  the  PI  controller  be  given  by 

Ullinit)  =  -KjX(t), 


where  Kf(t)  =  [kpf,  kif]T . 

The  resulting  closed  loop  system  is  the  desired  reference  system  for  adaptive  tracking 
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(37) 


(38) 


where  Am  is  Hurwitz.  The  total  control  input  for  subsystem  II  is  formed  as 

Uu(t)  =  umn(t)  +  ulad(t),  (39) 

where  Uu(t)  =  u\(t). 

Let  e(t)  =  X(t)  —  Xm(t)  be  the  tracking  error  signal.  The  tracking  error  dynamics  can  be  formed  as 

e(t)  =  Ame(t)  +  b(Akj  (t)X(t)),  e(0)  =  e0  =  0,  (40) 

where  A Kf(t)  =  Kf(t)  —  KJ  denote  the  parameter  errors.  Kj  are  ideal  values  of  the  adaptation  parameters. 
The  adaptive  control  can  be  designed  as  follows: 

ulad(t)  =  Kj(t)X(t),  (41) 
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with  the  following  adaptation  law: 


kf(t)  =  -YfX{t)eT(t)Pb,  Kf( 0)  =  Kf0,  (42) 

where  Tf  =  Tj  >  0  is  the  matrix  of  the  adaptation  gains. 

Defining  the  following  Lyapunov  function  candidate 

V{e(t),  A Kf(t))  =  eT(t)Pe(t )  +  A Kj (*),  rpA Kf(t),  (43) 

where  P  =  PT  >  0  solves  the  algebraic  Lyapunov  equation 

ATmP  +  PAm  =  -Q,  (44) 

for  arbitrary  Q  >  0.  The  time  derivative  of  the  Lyapunov  function  (43)  along  trajectories  (52),  (54)  is 

V(t)  =  -eT{t)Pe{t)  <  0.  (45) 

Hence  the  derivative  of  the  Lyapunov  function  is  negative  semidefinite,  therefore  all  signals  are  bounded. 
The  application  of  Barbalat’s  lemma  implies  asymptotic  convergence  of  tracking  error  to  zero. 

2.  Subsystem  I:  Fan  Control  Design 

The  goal  of  the  control  in  this  section  is  to  force  the  output  to  track  a  desired  command.  Output  is  the 
estimation  of  thrust  as  a  linear  function  of  the  plant  state.  The  controller  in  this  section  is  a  combination  of 
a  nominal  I  (Integral)  controller  and  an  adaptive  controller.  The  integral  controller  is  used  to  decrease  the 
steady  state  error,  and  also  to  generate  a  reference  model  for  adaptive  controller. 

We  design  a  full  state  feedback  adaptive  controller  so  that  y(t)  tracks  a  given  smooth  trajectory  yc(t), 
while  all  other  signals  remain  bounded. 

Here  we  introduce  the  integral  error  of  output  as  an  additional  state  into  the  system  dynamics  and  write 

y(t)  =  y{t)  -  yc(t )  =  c.6x(t)  +  d3.Su3(t)  -  yc(t).  (46) 

The  new  state  will  be  ^ 

yi(t)  =  [  y(r)dT.  (47) 

J  0 

Then  we  have 

Sx3(t)  =  y(t)  -  yc{t)  =  c.6xi(t)  +  d3.Su3(t)  -  yc(t),  &c3(0)  =  0,  (48) 

where  Sx3(t)  is  the  integration  of  output  thrust  tracking  error. 

The  total  control  input  for  subsystem  I  is  formed  as 

Ui(t)  =  uiin{t)  +  uad(t),  (49) 

where  Ui(t)  =  b2u2(t)  +  b3u3{t). 

Let  the  integral  controller  be  given  by 

Ulinit)  =  -kia5x3(t).  (50) 

Then  the  reference  model  will  be 

^3m(f)  =  C.8x\ m{t)  T  ^3 . {kia8x^  (t))  7/c(t),  &£3m(0)  =  0.  (51) 

Let  es(t)  =  Sxs(t)  —  5x% m(t)  be  the  tracking  error  signal.  The  tracking  error  dynamics  can  be  formed  as 

e3(t)  =  c.ei(t)  +  ds.kiae^(t)  +  d3.Aka(t)8x3(t),  e(0)  =  e0  =  0,  (52) 

where  A ka(t)  =  ka(t )  —  k*.  k*  is  the  ideal  value  of  the  adaptation  parameter. 

The  adaptive  control  can  be  designed  as  follows: 

uad(t)  =  ka(t)Sxs(t ),  (53) 
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with  the  following  adaptation  law: 


ka(t)  =  -~ia5x3(t)e3(t)sign(d3),  ka( 0)  =  kao,  (54) 

where  7a  >  0  is  the  adaptation  gain. 

Lyapunov  function  candidate  and  stability  analysis  for  this  controller  is  similar  to  the  previous  controller. 
The  stability  analysis  is  done  only  for  each  subsystem.  The  complete  stability  analysis  for  the  decentralized 
interconnected  system  is  not  done  here. 


B.  Core  Engine  Adaptive  Control  Results 

The  control  objective  is  to  design  a  robust  adaptive  control  such  that  the  state  of  the  nominal  system 
tracks  a  smooth  trajectory.  The  linearized  nominal  system  for  this  simulation  is  x(t)  =  a.5x(t )  +  b.5u\(t)  = 
—  1  .SSx(t)  +  1. 74605?/ 1  (t).  In  this  simulation  we  assumed  that  torque  is  a  function  of  time  and  defined  it  as 
T(t)  =  1  +  0.1  sin(1.3(l  +  sin(4t))t).  State  tracking  error  is  defined  as  e(t)  =  x(t)  —  xm(t). 

The  control  parameters  and  initial  conditions  are 


Kf= 


9 

75 

,  Tx= 

5  0 

0  30 

,  Q= 

2  0 

0  2 

,  p= 

6.5133 

-1.0000 

-1.0000 

0.2035 

,  Kf(  0)= 

-9 

-75 

(55) 


In  this  simulation  rc(t)  is  a  step  signal  which  defines  the  reference  spool  speed  to  change  from  1  to 
0.95.  Simulation  results  are  shown  in  figures  15  to  18.  Figure  15  shows  spool  speed  history  for  nominal  and 
reference  plant  and  also  state  tracking  error  history.  It  can  be  seen  that  state  follows  the  reference  trajectory 
very  closely.  Figure  16  shows  fuel  flow  history  as  control  input.  Figure  17  shows  the  history  of  the  two 
adaptation  parameters  used  in  the  adaptive  control  design.  Figure  18  shows  excess  power,  core  power,  SFC, 
turbine  temperature,  core  pressure  ratio,  and  core  airflow  time  histories. 


1.02 

1 


0.96 

0.94 


- Ref.  Signal 

- 

2  3  4 

time  (sec) 

Integration  of  state  tracking  error 


Fuel  flow  control  input 


time  (sec) 


Figure  15.  Spool  speed  and  state  tracking  error  histo¬ 
ries 


Figure  16.  Fuel  flow  control  input  history 


In  figure  6,  a  simulation  of  a  step  change  in  reference  speed  from  100  to  95  percent  is  shown  for  the 
engine  core  being  controlled  by  the  linear  PI  controller  used  to  define  the  adaptive  reference  model.  The 
controller  is  attempting  to  maintain  the  references  speed  under  the  influences  of  a  varying  torque  input. 
Figure  15  shows  the  same  simulation  for  the  adaptive  controller,  using  the  same  time  trace  of  disturbance 
torque.  In  both  cases,  the  controllers  follow  the  reference  rather  quickly,  but  the  adaptive  controller  exhibits 
a  smaller  amount  of  variation  due  to  the  disturbance  torque.  The  controller  is  able  to  adapt  to  changes  in  the 
linearization  caused  by  nonlinear  characteristics  of  the  engine  as  it  moves  away  from  the  design  condition. 
This  is  highlights  one  of  the  key  reasons  for  the  choice  of  an  adaptive  control  structure  in  the  distributed 
architecture  design. 
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Adaptation  parameters 
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Figure  18.  Turbine  temperature,  SFC,  pressure  ratio, 
Figure  17.  Adaptation  parameters  histories  power  and  airflow  histories 


C.  Turbofan  System  Decentralized  Adaptive  Control  Results 

Adaptive  decentralized  turbofan  engine  control  using  fuel  flow  and  fan  vane  angle  is  done  in  this  section. 
The  decentralized  structure  has  two  control  loops,  one  for  the  fan  subsystem  and  the  other  one  for  the  core 
engine.  For  fan  subsystem  we  let  the  fan  exit  area  to  be  constant  and  equal  to  its  design  value  (i.e.  v,2(t)  =  1); 
hence,  the  vane  angle  is  the  main  control  input.  We  measure  the  spool  speed  and  estimate  the  thrust  and 
use  these  two  values  to  design  an  adaptive  controller  which  actuates  fan  vanes  angle.  The  main  goal  for  this 
loop  is  to  force  the  thrust  as  the  output  to  track  a  desired  trajectory.  For  the  core  engine  subsystem  the 
input  is  fuel  flow. 

The  control  parameters  and  initial  values  for  subsystem  II  are 


Kf= 


9 

75 

r  — 

,  L  X  — 

5  0 

0  30 

,  Q= 

2  0 

0  2 

,  p= 

6.5133 

-1.0000 

-1.0000 

0.2035 

,  Kf(  0)= 

-9 

-75 

(56) 


The  plant  initial  conditions  are  x(0)  =  xm(0)  =  1.00.  The  control  parameters  and  initial  values  for 
subsystem  I  are 

kia  =  -10,  7a  =  10,  fca(0)  =  10.  (57) 

In  this  simulation  yc(t)  is  a  step  signal  which  defines  the  reference  thrust  to  change  from  1  to  0.95. 
Simulation  results  are  shown  in  figures  19  to  23.  Figure  19  shows  the  history  of  the  load  on  the  core  engine, 
and  also  actual,  estimated  and  reference, thrust  time  histories.  As  it  is  apparent  actual  and  estimated  thrust 
both  follow  reference  signal  closely.  Figure  20  shows  spool  speed  history  for  nominal  and  reference  plant,  and 
also  time  history  of  integration  of  the  output  tracking  error  (i.e.  yj).  State  follows  the  reference  trajectory 
closely.  Figure  21  shows  fuel  flow,  vane  angle  and  exit  area  histories  as  control  inputs.  Figure  22  shows  the 
history  of  the  adaptation  parameters  for  the  two  adaptive  controllers  for  subsystems  I  and  II.  The  adaptation 
parameters  converge  to  their  steady  state  values  fast  enough.  Figure  23  shows  excess  power,  core  power, 
SFC,  turbine  temperature,  core  pressure  ratio,  fan  pressure  ratio,  core  airflow  and  fan  airflow  time  histories. 

The  simulations  indicate  that  the  decentralized  adaptive  controller  appears  to  be  operating  properly, 
although  no  formal  stability  proof  is  currently  available. 

Figures  10  and  19  show  the  simulations  of  the  complete  distributed  control  system  using  both  the  linear 
controllers,  from  which  the  reference  models  are  based,  as  well  as  the  full  adaptive  controllers.  The  figure 
shows  a  step  change  in  requested  thrust  from  100  to  95  percent.  No  disturbance  inputs  have  been  added 
to  the  simulation.  The  simulation  shows  accurate  thrust  tracking  under  both  techniques,  with  the  adaptive 
control  excelling  in  overshoot  prevention.  As  seen  with  the  core  controller  alone,  the  adaptive  controller 
shows  increased  robustness  to  variations  in  the  model  caused  by  nonlinearities  away  from  the  design  point. 
This  robustness  to  model  variation  shows  promise  to  improving  performance  the  various  issues  that  may 
arise  in  a  non-ideal  distributed  environment,  such  as  communication  interruption  or  delay. 
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V.  Conclusion 

For  distributed  control  purpose  we  developed  a  simplified  dynamic  model  of  a  generic  turboshaft  engine 
driving  a  lift  fan.  Core  and  fan  subsystems  are  defined  as  the  subsystems  which  we  intend  to  control  using 
a  spatially  distributed  control  system.  In  this  distributed  model  fuel  flow  is  the  control  input  for  the  core 
subsystem,  and  vane  angle  is  the  control  input  for  the  fan  subsystem.  For  this  model  we  developed  a 
distributed  linear  control  structure.  Afterwards  we  used  the  linear  controllers  to  construct  desired  reference 
models  for  our  decentralized  adaptive  controller.  The  simulations  show  accurate  spool  speed  tracking  and 
thrust  tracking  under  both  linear  and  adaptive  techniques,  with  the  adaptive  control  excelling  in  overshoot 
prevention  with  less  oscillation  in  tracking.  The  current  competitors  to  Distributed  Adaptive  Engine  control 
are  the  traditional  FADEC  scheme  and  simple  distributed  control.  Both  Distributed  Adaptive  Engine  Control 
and  simple  distributed  control  offer  a  substantial  improvement  in  the  cost,  weight,  and  modularity  over  the 
traditional  FADEC  approach.  Distributed  Adaptive  control  has  additional  benefits  in  system  reliability, 
controller  stability,  and  modularity  of  the  entire  engine  system  that  give  it  a  significant  advantage  over 
any  simple  distributed  controller  currently  under  development.  In  future,  we  will  investigate  the  effect 
of  networked  control  system  issues  such  as  network  induced  delay  and  packet  drop  on  the  stability  and 
performance  of  the  developed  distributed  control  structure  for  turbofan  systems. 
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Overview 


Ieoh  Eco^Ei^^ifii^ 


•  Objective:  Research  the  control-theoretic  issues  associated  with 
distributed  control  of  non-homogenous  systems  containing  adaptive 
components.  Develop  the  functional  architecture  of  a  distributed 
control  system  in  a  gas  turbine  application. 

•  Approach:  Georgia  Tech  will  research  and  assess  approaches  to 
distributed  control  while  Aurora  identifies  applications  and  builds 
models.  The  control  schemes  developed  by  Georgia  Tech  will  be 
assessed  on  a  generic  gas  turbine  engine  model  coupled  with  a 
communication  emulator. 

•  Goals:  Demonstrate  the  stability  of  a  distributed  controller  in  a  Gas 
Turbine  and  identify  potential  performance  benefits  of  this  new  control 
architecture 
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Impact  of  Distributed  Engine  Control  on 
Future  Engines 


Current  Engine  Control  Architecture  is  Outdated 


I 

i 

Ieoh 


>  Essentially,  FADEC  is  a  fuel  controller  that  has  been  made  digital 

>  Complexity  has  increased  with  number  of  sensors/actuators,  need  for  health  management, 
etc. 

>  Further  advancements  are  hampered  by  unfavorable  complexity/reliability,  cost/benefit 
ratios 

>  A  change  in  overall  architecture  is  needed 


•  Smart  Sensors/Actuators  and  Smart  Components  Offer  a  New  Paradigm  for 
Future  Engines 

>  Dramatic  reduction  in  wiring  /  complexity 

>  Distribution  of  health  management  functions  to  smart  elements 

>  Component  performance  improvements 

>  Better  fault  handling 

>  Flexibility  to  incorporate  performance-improving  component  tailoring 

>  Modularity  of  control  system  and  engine  components 
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Turbofan  Model  Assumptions  and 
Simplifications 


Simplified  First-principles  based  dynamic  model  used  for 

control  law  development 


•  Only  the  dynamics  of  the  spool  are  considered. 

•  The  analysis  assumes  sea  level  static  operation  of  the  engine  system,  typical  of  a 
hovering  condition  for  a  VTOL  UAV. 

•  The  turbine  expands  the  exhaust  gases  of  the  core  perfectly  to  ambient  pressure, 
such  that  all  possible  work  is  extracted  from  the  flow.  The  exhaust  of  the  core  is 
assumed  to  have  no  impact  on  the  thrust  produced  by  the  fan. 

•  Non-ideal  efficiencies  are  assumed  for  the  fan,  compressor,  combustor,  and  turbine, 
but  all  other  components  are  assumed  to  operate  ideally. 

•  The  core  is  modeled  as  a  simple  gas  generator  with  no  secondary  or  bleed  flows  for 
turbine  cooling  or  other  uses. 

•  The  mass  and  temperature  of  the  fuel  is  ignored  in  the  calculations.  The  fuel  is  simply 
modeled  as  adding  heat  to  the  thermodynamic  cycle  in  accordance  with  its  mass  flow 
rate  and  heat  of  combustion. 

•  The  calculations  assume  that  the  specific  heat  of  the  air  going  through  both  the  core 
and  the  fan  is  independent  of  temperature. 
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f 
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Inputs: 

i:  fan  blade  incidence 
o:  nozzle  position 
q:  incremental  fuel  flow 


The  dynamics  of  the  spool  can  be  described  from  Newton's  second  law  by  the  sum 
of  the  torques  produced  by  the  turbomachinary  components  (T ),  as  well  as  the 
overall  rotational  inertia  of  the  system  (I). 

/w  =  Tt  -  %  -  Tj 


t,  c,  and  f  represent  the  turbine,  compressor,  and  fan  respectively. 

Torque  and  power  for  each  component  determined  analytically  from 
thermodynamic  and  mechanical  matching  between  components 
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x(t):  state,  non-dimensional  spool  speed 

y(t):  output,  thrust 

u^t):  fuel  flow  control  input 

u2(t):  fan  exit  area  control  input 

u3(t):  fan  vane  angle  control  input 


£(t)  =  /(z(*),ui{i),u2{*)^3(*)) 

y(t)  =  g(x(t'hu2(t),us(t)) 

-^dc-sj  ^  ^dts  X des  ==  1,  ^ties  ==  [1?  fi  0]  y<des  =  1 

<5t(£)  =  a.Sx(t)  4-  61  .<5if  1  (i)  4-  b2.Su2{t)  4-  b^.6u^(t) 

Sy(t)  =  Thrust  =  c.5x(t)  d2.5u2(t)  4-  cfe.i5-U3(t) 

a  =  —4.63,  61  =  1.78,  b-2  =  0.31,  63  =  2.96.  c  =  1,921  di  =  0,  d2  =  0.215,  ds  =  -1.44 
«i (t)  €  [0.3,  1.2]  u2(t)  G  [0.8,  1.2]  ti3(f)  €  [-10,  10 }(de.g) 
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Distributed  Control  Concepts  for  Turbofan 
Systems 
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Control  Approaches  for  Gas  Turbines 
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Distributed  Control  Structure 


Ref.:  D.  Culley  and  A.  Behbahani,  Communication  Needs  Assessment 
for  Distributed  Turbine  Engine  Control,  AIAA  2008-5281 
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The  highlighted  subsystems  are 
candidates  for  simplified,  distributed 
control  structure. 


In  our  model,  system  I  is  analogous  to 
the  “Inlet  Fan”,  and  system  II  is 
analogous  to  other  three  components 
(i.e.,  Compressor,  Combustor,  Turbine 
Nozzle). 
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Possible  sensor  and  actuator  locations 


Temperature  sensors 
sensors 


k- 


Vane  angle  actuation 
Exit  area  actuation 
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Example  Distributed  Control  Architecture 


"Ieoh 


Fuel  flow  loop  closed  by  speed 
measurement 


Vane  and  Area  loops  closed  by  thrust 
estimation  and  speed  measurement 
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Decentralized  Linear  Controller 
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Decentralized  Linear  control  structure  for 
turbofan  system 


Georgia 
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x(t):  state,  spool  speed 
y(t):  output,  thrust 
uM(t):  fuel  flow  control  input 
u,(t):  fan  vane  angle  control 
input 


Reference  speed:  rc(t)  =  [yc(t)  -  ydes  -  d3(u3  -  U3des) W  180)]/c  +  xde3 
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Decentralized  Linear  Control 


Ieoh 


System  I:  Fan  Control  Design 

«a(<)  =  kayi(t)  y[t)  =  y(i)  -  yc(i)  yi(t)  =  f  y(r)dT 

Jo 

ha  <  0 


Obtain  speed  reference  signal  r_c(t)  using  y_c(t): 

rjt)  =  [yc(0  -  Vdes  ~  <h{uz  ~  «3des)(7r/180)]/c  +  xde3 


System  II:  Core  Control  Design 


ui[i)  =  kpfX  +  kifXjft)  =  ^1 W  -  rc(0:  xi (*)  =  f  x(r)dr 
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Simulation  -  Core  Linear  Control 


T(t)  =  1  +  0.1  sin( '1.3(1  +  sin(4t))£) 
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Torque 


Georgia 

Ieoh 
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fuel  control  input  history 


Torque  Thrust 


Simulation  -  Decentralized  Linear  Control 
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Simulation  -  Decentralized  Linear  Control 
(Cont.) 
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MS  H  7  BCiBNOKB 
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Simulation  -  Decentralized  Linear  Control 
(Cont.) 
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Distributed  Adaptive  Controller 
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Decentralized  adaptive  control  structure  for 
turbofan  system 
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x(t):  state,  spool  speed 
y(t):  output,  thrust 
uM(t):  fuel  flow  control  input 
u,(t):  fan  vane  angle  control 
input 


Reference  speed:  rc(f)  =  [yc(t)  -  ydes  -  w3des)  (tt/180)]/c  4-  xdes 
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Decentralized  adaptive  control  using 
combined  fuel  and  vane  angle  control  inputs 


System  II:  Core  Control  Design 


Ieoh 

1 4w 


Plant  model: 


<5£i(t)  =  a.<foi(0  +  bi.&u\(t)r  i (0)  =  <5x2(0)  =  0. 

0  =  -  rc, 


Plant  model: 


Nominal  PI 
Control: 


Reference 

model: 


Sys  II  Controller 
(fuel  flow): 

Adaptation 

law: 


(5x2(0 


a  0 

1  0 
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<5x2  (f) 
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^  1  {^)  ~ l- 
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-1 


X(t) 


«liin(*)  =  ~KjX(t). 


-^/(O  —  [kpfi  ^i/]J 
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a  0 

"  h  ) 

a'/t) 
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^2m  (0 
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1  0 
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V 

(0 

-i 

rc{t). 


—s' 

Am 


U„(t)  =  uUi  „{t)  +  uuj(t)  ulad(t)  =  Kj  (t)X(t) 

K;(t)  =  -T/X(t)eT(t)Pb,  Kj  (0)  =  Kfo  r,  =  Vj  >  0 


N.  Hovakimyan,  Robust  Adaptive  Control  Course  Notes,  University  of  Illinois  at 
Urbana  Champaign  (UIUC),  2009. 
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Decentralized  adaptive  control  (Cont.) 
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System  I:  Fan  Engine  Control  Design 


Output  (Thrust):  =  C.5x{t)  +  d2.$U2(t)  4-  tfe-JlMO 

Output  Error:  y(t)  =  y(t)  -  yc(t)  =  c.Sx(i)  -  d2.Su2(t)  +  d^.5u^{t)  -  yc(t)  yj(t)  =  /  y(r)dr 

J{  i 


Plant  model:  6±s(t)  =  y(t)  -  yc(t)  =  c.5xi(t)  4-  d2.Su2(t)  -  yc(t ),  &r3(°)  =  0 


Reference 

model: 


^3m(0  =  c.<faim(t)  4-  d3.(M<tea(0)  -  yc(t ),  <fa3m(0)  =  0. 


Sys  I  Controller 
(vane  angle): 


Ui(t)  =  uun(t)  4-  tw(t)  kiaSx2(t)  ua(i(t)  —  ka(t)Sxs(t) 

ho.  <  o 


Adaptation 

law: 


MO  =  -la^(t)e3(t)sign(da),  ka(0)  =  kaG  7 a  >  0 


Ref.:  N.  Hovakimyan,  Robust  Adaptive  Control  Course  Notes,  University  of  Illinois 
at  Urbana  Champaign  (UIUC),  2009. 
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Simulation  -  Core  Adaptive  Control 
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Torque  Thrust 


Simulation  -  Decentralized  Adaptive  Control 
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Simulation  -  Decentralized  Adaptive  Control 
(Cont.) 
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Excess  Power 


Simulation  -  Decentralized  Adaptive  Control 
(Cont.) 
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Summary  and  Future  Work 
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•  Stated  Objective: 

>  Research  the  control-theoretic  issues  associated  with  distributed  control  of  non- 
homogenous  systems  containing  adaptive  components.  Develop  the  functional 
architecture  of  a  distributed  control  system  in  a  gas  turbine  application. 

•  Progress:  Demonstrated  distributed  adaptive  control  of  a 
simulated  turbofan  type  engine. 

>  Developed  dynamic  turbofan  engine  model 

>  Demonstrated  linear  distributed  control 

>  Demonstrated  adaptive  distributed  control  of  engine 

•  Future  goal  is  to  demonstrate  distributive  adaptive  control 
on  a  real  engine 

>  GE-80  developing  gas  turbine  driven  variable  pitch  fan 

•  Program  pending 
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•  Dealing  with  Networked  Control  System  (NCS)  issues 

>  Packets  dropout  during  transmission  and/or  incomplete  multiple- 
packet  transmission 

>  Network  induced  delays:  sensor-to-controller  delay  and 
controller-to-actuator  delay 

>  Bandwidth  limitations 

>  The  insertion  of  the  communication  network  in  the  feedback 
control  loop  makes  the  analysis  and  design  of  an  NCS  complex. 
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Future  Tasks  Required  for  Live  Engine  Test 
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Ieoh 
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•  Engine  control  laws  for  off-design  conditions 

>  Starting 

>  Steady-state  part-power  control  at  various  throttle  settings 

>  Full-scale  throttle  transients  (idle-to-max,  max-to-idle,  etc) 

>  Temperature,  pressure,  altitude  variation  accommodation 

•  Control  laws  for  fault  detection  and  safety 

>  Sensor/actuator  failure  recognition  and  accommodation 

>  Implementation  of  mitigation  strategies  for  over-speed,  surge,  etc 
to  prevent  catastrophic  hardware  damage 

•  Hardware  integration 

>  Controller/sensor/actuator  bandwidth  requirement  definition 

>  Real-time  embedded  system  details 
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Backup  Slides 


Noting  that  the  power  output  of  each  of  those  components  is  defined  by  V  =  Tw 

and  expressing  the  key  parameters  as  a  fraction  of  the  design  value,  previous  eqn  can  be  rewritten  as 


(  60  'Pdes  \  (  Vt  'Pudcs 

U*2  INlJ  \  Vt4es  Vdes 


Vr  V, 


c,des 


V 


c.des 


Vd 


S3 


t?-)  /  (£-) 

Pf,des/  \JVdes  ) 


where  rotational  speed  (N)  has  been  expressed  in  rpm  instead  of  rad/s. 

Compressor/fan  power,  written  as  a  fraction  of  the  design  power, 

j -  t-a 


W2 

^2,^63 


P-2  j  T2  VcMs  l  -  l 

p2,std  V  P‘2  ystd  *}c  ^  t  I 

1  ncHdes  ~  1 


Turbine  power  can  be  calculated  as  a  function  of  the  turbine  inlet  temperature.  Written  as  a 
fraction  of  the  design  power, 

Pt  wc  P2  j  T2  rjt  1  -  (x c7Tfc)  ^ 

Pi^des  ^^c,des  P2,std  V  P2  hstd  Vt^des  1  „  — .  ,  ) 

1  -1  \  n  c h de snb.de&  } 
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For  the  fan  and  compressor,  the  pressure  ratio  and  efficiency  of  each  component  is  expressed  as  a 
function  of  the  spool  speed  and  air  flow,  both  corrected  to  inlet  conditions.  The  change  in 
incidence  angle  of  the  fan  blades  from  the  design  condition  (delta-i)  is  also  input  to  the  fan  map. 


7TC  Tjc 
^c,des  Y}c,des 


=  Compressor  M  ap 
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^ des  V  'std 


rif 
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F  anM  ap 


For  the  purposes  of  the  Phase  I  investigation,  generic  maps  have  been  generated  for  the  fan  and 
compressor  based  on  assumed  phi/psi  characteristics  and  other  aspects  of  compressor  theory. 
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